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Single-file diffusion behaves as normal diffusion at small time and as anomalous subdiffusion 
at large time. These properties can be described by fractional Brownian motion with variable 
Hurst exponent or multifractional Brownian motion. We introduce a new stochastic process called 
Riemann-Liouville step fractional Brownian motion which can be regarded as a special case of 
multifractional Brownian motion with step function type of Hurst exponent tailored for single- 
file diffusion. Such a step fractional Brownian motion can be obtained as solution of fractional 
Langevin equation with zero damping. Various types of fractional Langevin equations and their 
generalizations are then considered to decide whether their solutions provide the correct description 
of the long and short time behaviors of single-file diffusion. The cases where dissipative memory 
kernel is a Dirac delta function, a power-law function, and a combination of both of these functions, 
are studied in detail. In addition to the case where the short time behavior of single-file diffusion 
behaves as normal diffusion, we also consider the possibility of the process that begins as ballistic 
motion. 
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I. INTRODUCTION 

Single-file diffusion (SFD) refers to the motion of par- 
ticles in quasi-one-dimensional channels and pores which 
are so narrow that the particles are unable to pass each 
other. The exclusion of a mutual passage of the diffus- 
ing particles means that the sequence of particle labels 
does not change over time. SFD is encountered in many 
physical, chemical and biological systems, which include 
the molecular and atomic motion in zeolites and nan- 
otubes, particle flows in microfluidic devices, ion trans- 
port in cell membranes, colloidal motion in narrow tubes, 

etc. p, a i is tag. 

The main feature of SFD is that for diffusion time 
t smaller than typical inter-particle collision time r c , 
the particles diffuse normally and satisfy Fick's law 
with its mean-square displacement (MSD) A 2 (t) := 

x(t) — x(0)] 2 ^) given by 



lim A 2 (t) 



2D t, 



(1) 



with Dq the diffusion coefficient. In other words, for 
t <C t c , the motion is just ordinary Brownian motion, 
which is a Markov process. However, for (>r c , 



lim A 2 {t) = 2FVi, 

t»T c 



(2) 
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where F is the SFD mobility. Recall that diffusion that 
does not satisfy Fick's law is known as anomalous diffu- 
sion with MSD satisfying A 2 (£) oc t a , a ^ 1. It is called 
superdiffusion when a > 1, and subdiffusion when a < 1. 
Thus, the long-time behavior of SFD belongs to anoma- 
lous subdiffusion, which is non-Markovian, indicating the 
motion is correlated. Note that SFD displays anomalous 
diffusion characteristics even when particle-channel in- 
teractions are not taken into account. 

Another way to characterize SFD is through its proba- 
bility density function (or propagator) P(x, t). The prob- 
ability of finding a particle at position x at time t, if it 
is initially at the origin, approaches the Gaussian propa- 
gator after a long time: 



P(x,t) 



1 



v/2vrA 2 (i) 



exp 



2A 2 (t) 



For (<t c , one has 



P(x,t) = 



and for t r c , 



1 



exp 



4£>nt 



P(x,t) = 



exp 



x 2 \ 



The notion of SFD was first introduced by Hodgkin 
and Keynes [8| who used it to describe the diffusion of 
ions through narrow channels in biological membranes. 
Harris was first to provide a theoretical derivation of (TTJ 
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and Q for SFD based on statistical argument Sub- 
sequently, this result was obtained using various models 
and methods by several authors including Levitt 1(311 . 
Fedders [ll|, van Beijeren et al [l2j and Karger 13j j. 
Recently, there are also attempts to model SFD based 
on fractional diffusion equations and fractional Langevin 
equations 

EE H E3 • Despite of the numerous the- 
oretical models and numerical simulations, experimental 
evidence for the occurrence of SFD was only obtained 
quite recently [H, El) H3, HH HI]. The main reason is 
that there is a lack of ideal experimentally accessible sin- 
gle file systems. 

The main aim of this paper is to propose some stochas- 
tic processes to describe the SFD. We do not address 
the detailed mechanism of SFD, but instead we empha- 
size more on the possibility of finding random processes 
which have the basic properties of SFD. Various types 
of fractional Langevin equations are considered in order 
to see whether they yield the stochastic processes which 
satisfy the basic statistical properties of SFD. 



by 



(B H (t)B H (s)) = (\t\ 2H + \s\ 2H -\t- s\ 2H ) , 



where 



C, 



H 



T(l - 2H) cos(vriJ) 



Bh (t) is not a stationary process, but its increment pro- 
cess is stationary. FBM is a self-similar process which 
satisfies for all a £ R+, 

B H (at)^a H B H (t), 

where = denotes equality in all finite distributions. The 
stationary property of the increments of B[{{t) allows the 
following harmonizable representation for the process: 



B H (t) = 



27T J-oo M 



II. MODELING SINGLE-FILE DIFFUSION BY 
STEP FRACTIONAL BROWNIAN MOTION 



In this section wc introduce a generalization of stan- 
dard fractional Brownian motion (FBM) [23] called step 
fractional Brownian motion (SFBM) and show that it 
can be used to describe the basic statistical properties of 
SFD. FBM has been widely used to model many areas 
such as turbulence, internet traffic, financial time series, 
biomedical processes, etc. One limitation of FBM model 
is that the long time (or low frequency) behavior that ex- 
hibit long-range dependence, and the short time (or high 
frequency) behavior that characterizes fractal property 
are both described by a single Hurst parameter H. Fur- 
thermore, a constant Hurst parameter is too restrictive 
for many applications. During the past decade different 
generalizations of FBM have been proposed to address 
this problem. Among them the most well-known is the 
multifractional Brownian motion (MBM), which was in- 
troduced independently in [24| and For MBM the 
Hurst parameter H is replaced by H(t), a deterministic 
function depending on time. MBM was later extended to 
generalized multifractional Brownian motion (GMBM) 
in order to model systems which require H(t) to be a 
irregular function of time [26|, [2?}. However, there exist 
processes and phenomena which exhibit abrupt changes 
of Hurst parameter requiring H{t) to be a piecewise con- 
stant function of time. For the description of such behav- 
ior, Benassi et al [H], Hi| introduced the step fractional 
Brownian motion (SFBM). A similar process known as 
multiscale fractional Brownian motion with its Hurst pa- 
rameter varying as a piecewise function of frequencies 
was also studied by several authors [H, [3l| . 

Recall that the standard FBM Bfj(t) is a Gaussian 
process with mean zero and correlation function given 



where < H < 1, t £ K and fj(u>) is the Fourier transform 
of rj{t) — the standard white noise defined by 

(ry(t)> = 0, (77(4)77(5)) =8{t-s). 

For modeling a process that evolves from time t = 0, 
instead of using the usual or standard FBM (which begins 
at time t — —00), it will be more appropriate to use an 
alternative FBM that starts at time zero. This second 
type of FBM is known as Riemann-Liouville FBM (RL- 
FBM), which is defined as the RL fractional integral of 
white noise 1321: 



w H {t) 



1 



T(H+\) Jo 



(t-u) H ^^r](u)du, (eK + , H > 0. 

(3) 



Wjy(i) is a Gaussian process with zero mean (Wjr(i)) = 
and correlation function given by 



C WH (t,s) = (W H (t)W H (s)) 

t H-l/2 s H+l/2 _ / 1 

(H+l)T(H+^ ' * V2 



H,1,H + 



2' t 



when s < t. Here 2-Fi(a, 6, c; z) denotes the Gauss hyper- 
geometric function. The variance of the process Wn{t) 
is 



Var(W H (t)) - (W H (t) 2 ) = 



2H 



2HT (H+\Y 



(4) 



Note that for the standard FBM B H (t), the Hurst 
parameter H should lie in the range (0, 1), whereas for 
RL-FBM Wn{t), H takes any positive real value. Both 
Bn{t) and Wu(i) reduce to ordinary Brownian motion 
when H = 1/2. In contrast to Bfj{t) which has sta- 
tionary increments, the increments of Wn(t) are non- 
stationary. Due to the failure of its increments to be 
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stationary, Wh (t) does not has a harmonizable represen- 
tation. When t — ► oo, RL-FBM approaches the standard 
FBM [32|. 

Now we want to consider the step fractional Brown- 
ian motion (SFBM). Such a generalization of FBM was 
first introduced for standard FBM using the harmoniz- 
able representation as follow [28|: 



by Hi > and H2 > 0. For this simple case, we can 
write 



l[o,r)(t) 

r(#i + ±)7o 

1[t,oo)(*) 



(t - u) Hl ~^ri(u)du 

(8) 

(t - u) H2 ^^r](u)du. 



1 



V^J- 



where 



(5) 



r (if 2 + 1) Jo 

From the above definition, one sees that Wh u h 3 (£) is a 
Gaussian process with zero mean and correlation function 

(W Hl ,H a {t)W HltHa (8)) 

=l[o,T)(*)l[o,T)(s)Ci?i,ifi(* ; s) 

+ l[r,oo)(i)l[0,T)( s )CHi,H 2 (* 5 s ) 

+ l[r,Oo)(*)l[T,CX))( S )Cj?2,J?2(*) s) 



with To = —00 and tjv = 00, Hi £ (0, 1), 1/(4) = 1 if t € I 
and = if t £ I. This is an adaptation of MBM, 

which is defined for a time-dependent Hurst parameter 
H(t). Due to the absence of such representation in RL- 
FBM, we generalize W H {t) to RL-SFBM based on the 
moving average representation ([3]): 



W H[t) (t) 



1 



r(H(t) 



when s < t, and Ga^Hi (*i s )j = 1, 2, is given by (J7]) 
with H(t) = Hi and i?(s) = Hj. Similarly, its variance 
(W Hl ,H 2 (t) 2 ) is given by © with i? replaced by if(t). 
Thus, Wh 1 ,h 2 the RL-SFBM with two scales, behaves 
like Wn-iif) for i G [0, r), and behaves like W/r 2 (i) when 
t £ [t, 00) (see FIG. [1]). W^i^t) is piecewise self- 
similar, it is self-similar of order Hi in the time interval 
I (t- u) H ^^~ 1 / 2 ?7(u)(i-u, (6) [0, r), and of order iJ 2 in the time interval [r, 00). 

In order to use RL-SFBM for modeling SFD, it is nec- 
essary to carry out some minor modifications to the def- 
inition of WH 1 .H 2 {t)- © an d © are to be replaced by 



2) J o 



with H(t) the piecewise function given by except 
that in this case To = 0, and Hi £ (0, 00). Its covariance 
is given by 



C WH( . ) (t,8) = (W m .)(t)W Hi . ) (a)) 

t H(t)-l/2 s H(s)+l/2 

= r(H( S ) + l)T(H(t) + ±) 
x 2 Fi(±-H(t),l,H(s) + ^^ 



W H (t)(t) 



XH(t) 



r (H(t) + i) J 



(t - u) H ^-^ 2 r)(u)du, (9) 



and 



(7) 



iy gl|g| (t)= g^ [\t-u) H ^ v { u )du 



(t - u) H2 ~^ri{u)du, 



r (#1 + 1) 7o 

XH 2 1[t,oo)(*) /" 



if s < i, and its variance (W^.) (i) 2 ) is given by (|4]) with 
-ff replaced by 

The property of global self-similarity does not apply to 
both MBM and SFBM. In the case of MBM, the notion 
of self-similarity is replaced by the local asymptotic self- 
similarity L 2j|, which is also satisfied by SFBM and RL- 
SFBM with some modification. Suppose W H ^{t) is a 
RL-SFBM with scaling function Hit) defined above. For 
all t £ (Ti_i,Tj), 



r(ff 2 + i) Jo 



(10) 



lim 

e^0 



W H {.){t + eu)-W H (.){t) 



The convergence is in the sense of distributions. In other 
words, the tangent process of RL-SFBM for each scale 
Hi is B Hi , a FBM indexed by Hi. 

For modeling SFD, we use RL-SFBM with single 
change of scale, that is the process ((6J with H(t) — 
Hilm T \ + i?2l[ r ,oo)- To be more specific, we denote this 
process by Ws 1 ,H 2 {t), the two-scale RL-SFBM indexed 



where XHn i = 1)2, are positive constants which are 
introduced for the purpose of obtaining the correct coef- 
ficients for the MSD of the diffusing particles. Note that 
by definition, {W H (■)($) = and W H( . } (0) = 0. There- 
fore the variance of Wjj^(t) is equal to the MSD. 

Now we want to see how the RL-SFBM with two scales 
can be used to described the basic properties of SFD. By 
letting Hi = 1/2, H2 = 1/4, and XHi in terms of diffusion 
coefficient Do and SFD mobility F with xh ± = 2Do, 
X h 2 = r(0.75) 2 F, then the MSD of W Hl ,H 2 (t) d) is 
equal to 



(w HuH2 (t) 2 ) = ([W Hl H 2 {t) - W Hl H 2 (0)Y 



i2D t, 
\2FV~t, 



t£[0,r) 
t £ [t,oo), 



which is the same as the MSD of SFD. The value of r is 
set to equal to y/F/Do so that the MSD is continuous. 
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FIG. 1: Two step RL-SFBM W Hl ,H 2 (t) with Hi = 0.5, 
H 2 = 0.25, t = 5, xh x = 1 and Xh 2 = \/E>r(0.75) 2 /2. The 
smaller window shows the MSD of the process. 




t 



FIG. 2: Piecewise linear RL-MBM W H{t) (i) with 
H(t) = 0.5, for t £ [0, 3], H(t) = 0.25 for t £ [6, 10] and H(t) 
is linear in for t £ [3, 6]. XhU) — 1 f° r t £ [0, 3] and 
X H(t) = 2 for t 6 [6, 10]. For i £ [3, 6], X ff(t) is linear. The 
smaller window shows the MSD of the process. 



A sample path of the process Wh 1 h 2 (t) is shown in FIG. 

m 

There are some further generalizations of RL-SFBM 
which can be used to model SFD. For example, one can 
use a piecewise linear function H(t) in ^ so that H{t) = 
1/2 for t £ [0,ti), H(t) = 1/4 for t £ [r 2 ,oo) and H(t) 
is a linear function interpolating the points (n, 1/2) and 
(r 2 ,l/4) in the interval [ti,t 2 ]. If Xff(t) = 2-Do for t £ 
[0,n), Xff(t) = r(0.75) 2 F for t £ [r 2 ,oo), and XH(t) is a 
linear function of t in the interval [ti,t 2 ] so that XHU) 
is continuous, then the process WhU) (t) is a special case 
of RL-MBM, which provides a model for SFD that has 
continuous sample paths (see FIG. [5]). 

As remarked by Karger [2|, random walk model can 
only be regarded as an approximation to the real SFD 
systems. At short times, such systems first undergo bal- 
listic motion with the MSD A(<) ~ t 2 . In other words, 



0.5 
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FIG. 3: Three step RL-SFBM W Hl ,H 2 ,H 3 (t) with H x = 1, 

H 2 = 0.5, H 3 = 0.25, n =1,t 2 = 5, X h ± = 1, 
X 2 H2 = l/(2r(1.5) 2 ) and xl, = \/5r(0.75) 2 /(4r(1.5) 2 ). The 
smaller window shows the MSD of the process. 



there is the possibility of the direct transition from the 
ballistic regime to the single-file regime. Such a tendency 
becomes more prominent with increasing concentration; 
and it has been demonstrated by molecular dynamical 
simulations [33|. For modeling SFD that is ballistic at 
small t, we can use RL-SFBM with three scales. Here we 
consider the case of a SFD process with three regimes: 
initial ballistic regime followed by the normal diffusion, 
and finally the single-file diffusion region. For such a pro- 
cess, we let the time-dependent Hurst exponent in © to 
be 

H(t) = Hil[ 0)T1 )(t) + H 2 l[ TltT2 )(t) + i7 3 l[ T2 00 )(t), 

where Hi = I, iJ 2 = 1/2 and H% = 1/4. An example of 
such a process is shown in FIG. [3] 



III. MODELING SINGLE-FILE DIFFUSION BY 
FRACTIONAL LANGEVIN EQUATIONS 

In this section we want to examine whether it is possi- 
ble to describe the basic characteristics of SFD based on 
the various types of fractional Langevin equations. We 
shall first show that how RL-SFBM can be obtained as 
the solution to the "free fractional Langevin equation" 
(fractional Langevin equation without damping). This 
will be followed by discussion on the fractional general- 
ized Langevin equation and its extensions. 



A. Fractional Langevin equation 

First note that the definition of RL-FBM ([3]) can be 
written as: 

W H {t) = /f +1/2 V (t) 7 (11) 
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where the Riemann-Liouville (RL) fractional integral qI" 
is given by [H, [H, [H : 



(oi?/)(*) 



1 



(t-u) Q_1 /(w)du. 



The Riemann-Liouville fractional derivative is de- 
fined as 



D? := D? Q I?- a (12) 
for n — 1 < a < n. In view of the property that 

i? D?f(t) = i?D? ir a f(t) 

t a-k 

k - r(a-* + l) 



=/(*) - £ ■ 



onr k f] (o), 



we can consider (fTTj) as the solution of the "free" frac- 
tional Langevin equation 

Q D? +1/2 w H (t) = nQt), 

subject to the initial condition ( o-D^ _1 ^ 2 H# )(0) = 0. In 
a similar way, we see that can be re-expressed as 



W Hl ,H 2 (t) = Ol, 



_ T H(t)+l/2 



1[0,t)(*) qI, 



Hi + 1/2 



1[t,oo)(0 qI, 



H-2 + 1/2 



r)(t), 



and Wffi.HaC*) can b e regarded as the solution of 



D? W+1/2 W H{t) (t)= v (t), 



with jEf(t) = i?il[ . T ) + i?2l[r,oo)) an( l subject to the 
initial condition (oDf {t) ~ l,2 W H {t)\ (0) = 0. 

Recall that Brownian motion can also be regarded as 
the position process associated with the solution (veloc- 
ity process) of the usual Langevin equation. It is nat- 
ural to ask if a similar link exist for RL-FBM and RL- 
SFBM with fractional Langevin equation. For this pur- 
pose we first consider the following general type of frac- 
tional Langevin equation with two different fractional or- 
ders a and 7 [3^ |: 

(oflf+A) 1 ^^^^), 0<a<l,7>0, (13) 

where A > is the dissipative parameter, and r)(t) is the 
standard white noise. Here we remark that (qD" + A) 7 
can be regarded as "shifted" fractional derivative as com- 
pared with the un-shifted one qD" 1 . By using binomial 
expansion, the shifted fractional derivative can be for- 
mally expressed as 



00 
3=0 V ' 



7 U J oDt o(7_i) - 



Special case of (TT51) w ith a > 0,7 = 1 has been consid- 
ered previously in [39l l4p|. Ell. |42| : and the solutions of the 



fractional Langevin equation with a = 1, 7 > have also 
been studied in [4^, |4J] . It is found that for the general 
case, if one considers the solution v an (t) to (fT3"| as the 
velocity process, then the corresponding position process 
x an {t) depends only on the differential relation between 
v a .j(t) and x a .j(t) [45| . If the usual velocity-position re- 
lation is used, that is velocity is the ordinary derivative 
of position, then the variance of the position process does 
not depend on the long-time behavior of the correlation 
of the velocity process. One always get the variance of 
x an {t) behaves like Var(x an (t)) ~ t, just like the case 
of normal diffusion. The long-time dependence of the 
covariance of v an (t) which varies as i""" 1 does not en- 
ter in the leading term of the variance of x Q . 7 (t), it only 
appears as second leading term 45:]. 
On the other hand, if one assumes 



y(t) = oDtX an] p(t), fo-Df 1 Xa, 1 ;pj (0) 



0. 



where < j3 < 1, such that 



then we have shown in |45[ that Var(x a , Ti p(Jb)) ~ i?^~ x . 
Therefore for j3 — 3/4, one gets the correct long time be- 
havior for the MSD of SFD. Here we have used the fact 
that £c<,7;/3(0)=0. As for the short time behavior, it can 
be shown that Var(x a , Ti/3 (t)) ~ t 2 «7+2/3-i g|_ There- 
fore, in order to obtain the correct short time behavior 
for the MSD of SFD, we require cry = 1 - j3 = 1/4. We 
notice that for the characterization of the mean square 
displacements of SFD, the parameters a and 7 appears in 
the combination cry. Therefore, we can restrict ourselves 
to the case 7 = 1, for then the process x a ^.p(t) satisfies 
the following fractional Langevin equation: 



oDt x ot,i;p{t) =v a ,i{t) 
D?v ail (t) + Xv ail (t) = 77(f), 



(14) 



with initial conditions (qD^: x a -i-p)(0) = and 
(o-Dfu a; i)(0) = 0. Setting a = 1/4 and (3 = 3/4 re- 
covers the basic properties for SFD. On the other hand, 
setting a — j3 = 3/4 gives a process x a ,i;/3 that is ballistic 
(i.e., its MSD is ~ t 2 ) at small t, and sub-diffusive with 
exponent 1/2 when t is large. 

Notice that for n — 1 < a < n, the Laplace transform of 
the Riemann-Liouville fractional derivative o-D" is given 
by HI: 

n-l 

^Dff(s) = s a f(s) -J2* k [oDf-"- 1 /] (0), (15) 

fc=0 

where f(s) denotes Laplace transform of f(t). In solving 
the fractional Langevin equation (|13[) . we have assumed 
that (o-D" _1 Wq i7 )(0) = so that the Laplace transform 



of the solution u ai7 (s) satisfies 



fj(s) 



(s a + A) 7 ' 
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FIG. 4: A sample path of Xa^pit) where a = 1/4 and 
j3 = 3/4. The smaller window shows the MSD of the process. 



when n — 1 < a < n. The initial conditions that need 
to be specified now are the values of the ordinary deriva- 
tives of / at t = 0, which have natural physical inter- 
pretations. For the system (fT4"]) . we have assumed that 
( o D%x a;1; p)(0) = and ( D?v a; x)(0) = 0, which give the 
relation 

S^XaS-pis) = V a ,l(s), 
S a V a ,l(s) + Xv at i(s) = fj(s) 

for the Laplace transforms of x a ,i-,/3(t) and v at \(t). Com- 
paring the Laplace transforms of the Riemann-Liouville 
fractional derivative and Caputo fractional derivative 
(IT5l) and (fl6|) . we find that we can also interpret our 
solution to (|14[) as satisfying 



D?v a<1 (t) + \v a<1 {t)=r)(t), 



(17) 




FIG. 5: A sample path of x a ,i-,p(t) where a = 3/4 and 
13 = 3/4. The smaller window shows the MSD of the process. 



and the solution w Q . 7 (i) is obtained by taking the inverse 
Laplace transform of this equation. From the practi- 
cal point of view, the applicability of Riemann-Liouville 
fractional derivative is limited by the absence of phys- 
ical interpretation of the initial condition of the type 



r )(0) = vq, when a is not an integer. There is 



another definition of fractional derivative called Caputo 
fractional derivative which is defined as 



when n — 1 < a < n. The difference between the 
Caputo fractional derivative and the Riemann-Liouville 
fractional derivative p2[) lies in the order of taking dif- 
ferentiation and integration. In contrast to (|15[) . the 
Laplace transform of Caputo fractional derivative is given 
by HI 



$D?f{s) = 3 a f{8) 



k=0 



WO), (16) 



with initial conditions x a ^fj{Q) — and w Qi i(0) = 0, 
where now Caputo fractional derivative is used. In the 
following, we are only going to use Caputo fractional 
derivative. Therefore, we are going to use the symbol 
o-D" instead of the symbol q Df for Caputo fractional 
derivative and it should not incur any confusion. The 
short time behavior of MSD is sensitive to the initial con- 
ditions. For example, the stochastic process satisfying 
(fTTf with a = j3 = 1 and initial conditions v(0) ^ has 
MSD behaves like - t 2 as t -> 0. Whereas for v(0) = 0, 
the MSD behaves like ~ t 3 as t — > 0. We will consider 
the more general case where a;(0) ^ and v(0) ^ in 
the subsequent discussions. 



B. Fractional Generalized Langevin Equation 

We consider the generalized Langevin equation: 
D t x(t) = v(t), 



D t v(t) + 7(t - u)v{u)du = F(t), 



(18) 



where j(t) is the dissipative memory kernel, and F(t) is 
a Gaussian random force with zero mean and correlation 



(F(t)F(s))=C F (\t-s\). 



(19) 



When 7 (t) = XS(t) and C F (\t\) = S(t), d]) reduces 
to the ordinary Langevin equation. It has been shown 
that it is possible to obtain the position process x{t) as 
an anomalous diffusion if F(t) is considered as internal 
noise with long-tailed correlation (46|. More precisely, 
if 7(i) = \t~ K , where < k < 1, and fluctuation- 
dissipation theorem holds with Cp{t) = kBT^/lt), one 
then has ([x(t) — a;(0)] 2 ) ~ t K as t — > oo. On the other 
hand, if F(t) is an external noise with Cp(t) = cg\t\~ e , 
1< 6 < 1, one has (x(t) 2 ) - t 2K - g as t -> oo if 2k > 6 
(13, IHl . We thus have subdiffusion when < 2k— 6 < 1, 
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superdiffusion when 2k — 9 > 1, and normal diffusion for 
2k — 9 = 1. However, the short time behavior of x(t) is 
always ballistic (IrJ ]. 

In this subsection we want to examine whether it 
is possible to describe the basic characteristics of SFD 
based on the fractional Langevin equation in the follow- 
ing general setting: 

Dfx(t)=v(t), O<0<1, 



D?v(t) + j(t- u)v(u)du = F(t), < a < 1, 
Jo 

(20) 

where now ^Df and o-Df are Caputo fractional deriva- 
tives, 7(i) is the frictional kernel and F(t) is a Gaussian 
noise with zero mean and the following correlation: 

(F(t)F(s)) =C F (t- s ) = c e \t-s\- e . 

When (3 = 1, we have ordinary velocity, and for /3 ^ 1, 
the velocity is a fractional velocity. 

If the dissipative memory kernel is given by ^(t) = 
\t~ K , < K < 1, the second equation in (|20|) can be 
written in a more compact form: 



D?v{t)+ Xo liv{t)=F{t) 



(21) 



where ( = 1 — n and x — -T(C)A. Formally, when £ 
and F(i) = r](t), (|2T|) reduces to ordinary Langevin 
equation. Laplace transforms of ([20| and (|2Tj) give 



s^a;(s) — s' 3 1 io = £(s), 



s a v(s) - s a - L v + Xs~ C v(s) = F(s), 



(22) 



where x = x(0) and vq — v(0). We assume that v n ^ 0. 
From (f2"2"| one gets 



i/(s) 



x(s) = — 



F(s) v 



s s Q +^ (1 + x s ~ C ~ a ) sl3+1 (1 + X S ~ <: ~ Q )' 

(23) 

From now on, we only concentrate on the solution to x(t). 
Inverse Laplace transform of (f2"3"|) gives 

i(t) - x + votPE a+Ctf i +1 {~xt a+C ) 

+ I (t- ur+^Ea+^+p (- X (t - u) a+ <) F{u)du, 



E^(z)=J2 



where 



k=0 



T{nk + v) 



is the two-parameter Mittag-Leffler function [49( . Taking 
expectation value using (|19[) . we have 



(x(t)) =x + vot?E a+Ctfl+1 (~ X t a+C ) 



and 



a 2 xx (t) = ([x(t) - (x(t))} 2 ) = 2c e r {1-6) f u 



2ct+2/3-0-l 



x E a+i>a+/3 (-xw Q+C ) ^a+c,a+/3-e+i (-X«" +C ) du, 

where we have used the following identities [35] : for v > 
and 6 < 1, 



(i - u)~V'- 1 .E M , 1 , (-xO<fe 
=r(i - 9)t»~ e E^ v _ e+1 . 
The MSD is given by 

A 2 W =aL(t) + ((^))-x ) 2 

Now by usin g th e asymptotic properties of Mittag- 
Lcmcr function 1491: 



E^(-z)=J2 



(-1) 



A:=l 
1 



r {v — /ifc) 



CO, 



^(-^=^+0(2), z^O, 
one gets for t — > 0, 

«£.(*) ~t 2 "+ 2 ^, 

'i 2 * 3 , if 2a > 



A 2 W 



with the assumption that 2a + 2/3 > 9 so that <r xx (t) has 
a finite limit as £ — * 0. 

For large time asymptotic behaviors, we have generi- 
cally 



■ t 20-2C-e^ if2/3-2C>0, 
lni, if 2/3- 2C = 0, 

constant, if 2/3 - 2( < 9. 



Therefore, for the MSD, if 2a > 9, the large time asymp- 
totic of A 2 (t) is the same as <J xx (t). However, if 2a < 9, 
then as t — > 00, 



A 2 W 



r t 2/3-2C-2 Qj if 2/3- 2C > 2a, 



1 constant, 



if 2/3 - 2C < 2a. 



Now consider the general case with < ( < 1. 
Fluctuation-dissipation theorem requires 



(F(t)F(s)) = c e \t - si 6 = k B T x 



no ' 



which gives = 1 — £. In order the MSD satisfies the 
properties of SFD, we require its asymptotic behavior 
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~ t when t — ► and ~ \ft when i — > oo. This gives two 
possibilities: 



Case I 2a > 1 



■C, /?=-, 2/5- C 



Case II 2a < 1 - C, 2a + 2y3 + £ - 1 = 1, 
2/3-2C-2a=i. 

Case I implies that £ = —1/2, which is a contradiction 
to C G [0, 1). For Case II, we find that 



_ 3 3C 
a ~ 8 ~ T' 



H 4 8 



(24) 



The conditions a € (0, 1], /3 € (0, 1] and 2a < 1 - ( imply 
that C e [0,1/2). In other words, for any ( G [0,1/2), 
define a and by ([24| . Then the process x(t) gives a 
correct description of SFD. 

If we assume that x(t) is ballistic instead of normal 
diffusive at small t, then the possibilities are 



Case I 2a > 1 — Cj = 1, 2/3 — £ — 1 



1 

2' 



Case II 2a < 1 - C, 2a + 2/3 + £ - 1 

1 

2' 



2/3- 2C- 2a 



Case I gives /3 
we find that 



1, C = 1/2 and a > 1/4. For Case II, 



3C 
4 ' 



0- 



4 8 



The condition <1 leads to £ < 1/2. However, the con- 
dition 2a < 1 — C gives £ > 1/2, which is a contradiction. 
Therefore, in order that (f2"0"|) gives a suitable model for 
x(t) which is ballistic at small t and sub-diffusive of ex- 
ponent 1/2 at large t, we need to set /3 = 1, £ = 1/2 and 
a can be any number between 1/4 and 1. When a = 1, 
this reduces to the model (|18p in the beginning of this 
section. 

Finally, let us remark on the case where vq = 0. In this 
case, A(i) 2 = cr xx (t). We then find that the properties 
of SFD are satisfied if 



1 

4 



c, 



for C G [0,1/4). When C = 0, this gives a = 1/ 4 and 
/3 = 3/4, which agrees with the result in Section Ull Al 
On the other hand, if we require the process to be ballistic 
at small t, then 



3 
4 



c, 







2 4 



for £ £ [0,3/4). C — gives a — — 3/4, which again 
agrees with the result in Section fill Al 



C. Extended Fractional Generalized Langevin 
Equation 

In this section, we stretch our Langevin approach to 
an even more general setting which includes the cases 
discussed in Section ITlI Al and Section lTlI Bl Consider the 
following extended version of the fractional generalized 
Langevin equation: 



D?x(t)=v(t), 
Q D?v(t) + Xv{t) + 



X 

no Jo 



(t - u) c ^ 1 v(u)du = F(t), 



(25) 



where now the dissipative memory kernel "f(t) is given by 



7(f) = 2XS(t) 



X 

r(C)' 



(26) 



with < C < 1- When a = (3 = 1 and C = 1/2, this 
system was proposed as a model to describe SFD for 
which the MSD is ballistic (i.e. ~ t 2 ) when t is small. 
When A = or x = 0, the equations (|25|) reduce to 
the equations (|2"D|) considered in Section Ull Bl Therefore 
here we assume that A ^ and x ^ 0. As in Section III, 
taking Laplace transforms give 



X[S> s s Q +^(l + As- Q +xs" c " a ) 



(27) 



sP +1 (1 + As-" + xs~ c ~ a ) ' 



In the following, we assume that vo ^ 0. Let Ki(t) and 
K2(t) be respectively the inverse Laplace transforms of 



1 



and 



1 



s a+/3 (1 + X s -a + X s-C-a) S 0+l (1 + As~ a + X^ C ^ a ) ' 

Then the inverse Laplace transform of (|27)) gives 

x(t) = x + [ K 1 (t-u)F(u)du + v K 2 (t). (28) 
Jo 

From this, we find that the variance of x{t) and the MSD 
is given respectively by 

o% x (t) =2 f [ K 1 {v)C F (u~v)K 1 (u)dvdu 
Jo Jo 

A 2 (t)=2 / K 1 (v)C F (u-v)K 1 (u)dvdu + v 2 K 2 (t) 2 . 
Jo Jo 

(29) 

The asymptotic behaviors of the functions Kx(t) and 
K 2 (t) at small and large t are studied in Appendix [X] 
The result is: as t — > 0, 



As t — > oo, 



K 2 (t) ~ IP. 
K 2 {t) ~ t^-C" 



(30) 
(31) 
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As in Section fill Bl we assume the following generalized 
fluctuation-dissipation theorem [50j : 



Cj (/) =l L1 T (2XS(t) + ^t^- 1 



Then the MSD J29]) can be rewritten as 



(32) 



A z (t) = v^K 2 (ty + 2k B TX / Kx{uydu 

Jo 

+ 2k B T-^- f f Kxiv^u-vf^Kxi^dvdu, 
1 (.tj Jo Jo 



(33) 



30)) then implies that its small time asymptotic behavior 



A 2 W 



t 2(3 7 if a > 1/2, 

t*»+2/»-i, if a < 1/2. 



For the large time asymptotic behavior of the MSD, (|3"Tj) 
and (|3"3"j) give 



A 2 W 



hit, 

constant, 
t constant, 



if 2a + C > 1 and 2(3 - ( > 1 
if 2a + C > 1 and 2/3 - C = 1 
if 2a + C > 1 and 2/3 - ( < 1 . 
if 2a + C < 1 and /3 - ( > a 
if 2a + ( < 1 and /3 - ( < a 

(34) 



The short time asymptotic behavior of the MSD is gov- 
erned by the term d(t) in the dissipative memory kernel 
and the long time asymptotic behavior is governed by 
the term tS~ l . This gives the general form of dissipative 
memory kernel (|26|) the advantage of being able to inter- 
polate between the two particular cases with A = and 
X = considered in the previous subsections. Now to 
satisfy the characteristics of SFD, there are a few possi- 
bilities: 

Case I a> i = ± 2/3 - C - 1 = \, 
Case II a < i, 2a + 2/3 - 1 = 1, 2a + C > 1, 



2/3-C-l 
1 



2' 



Case III a < -, 2a + 2/3 - 1 = 1, 2a + C < 1, 
2/3-2C~2a= i. 



The condition 2a + ( < 1 is then automatically satisfied. 
The conditions < a < 1/2 and < (3 < 1 then imply 
that < C < 3/4. In other words, for any C 6 [0,3/4), 
define a and /3 by (|35[) . Then we obtain a solution cc(t) 
to ([25]) which has the characteristics of SFD. Compared 
to the solution in Section fill Bl which only allow ( to lie 
in the range [0, 1/2), we find that the extended fractional 
generalized Langevin equation (|25p can be used to de- 
scribe SFD in a larger range of £. In both cases, the 
maximum value of a is 3/8. 

Next we consider the case addressed in the paper [5C| , 
where the MSD is ballistic for small t, and become sub- 
diffusive ~ ^/t when t is large enough. For this to happen, 
there are a few possibilities: 

Case I a > i, /3 = 1, 2/3 - £ - 1 = ^, 

Case II a < i, 2a + 2(3 - 1 = 2, 2a + ( > 1, 

2/J-C-l = |, 

Case III a < i, 2a + 2(3 - 1 = 2, 2a + ( < 1, 

2/3-2C-2a= i 

For the first case, we find that /3 = l,C=l/2 and the 
only restriction on a is a > 1/2. When a = 1, this is 
the case considered in Case II and Case III imply 

that a = 3/2-/3 > 1/2 since (3 < 1. But this violates 
the condition a < 1/2. In summary, to characterize the 
behavior of a system which is ballistic at small t and sub- 
diffusive at large t, one can use the extended generalized 
fractional Langevin model (|25[) with (3—1 (so that the 
velocity is normal derivative of position), £ = 1/2 and a 
any values between 1/2 and 1. Compare to the result of 
previous section, we find that when A = 0, the conditions 
C = 1/2 and (3=1 are the same, but the range of a 
is from 1/4 to 1, which has larger range of values as 
compared to the case when A ^ 0. 

Finally we remark on the case where vq = 0. In this 
case, A 2 (t) = cr^ x (t). Therefore, the small and large time 
asymptotic behaviors of the MSD for SFD is satisfied if 



a 



4 2' P 2 4' 
and the MSD is ballistic at small t if 



< C < 



1 



3 C 
a = . 

4 2' 



o < c< i. 



We have used the fact that a > 1/2 implies 2a > 1 > (. 
Case I gives ( = —1/2, which is a contradiction. The 
two equalities in Case II imply that 2a + ( = 1/2, which 
contradicts 2a + ( > 1. For Case III, the two equalities 
lead to 



3 _ c 

8 2' 



0- 



8 2 



(35) 



D. Further Generalizations 

The results above can be easily generalized to frac- 
tional Langevin equation with more general (nonlocal) 
dissipative memory kernel 7(f) which behaves like ~ 
t K_1 , k € [0,1) when t — > and behaves like ~ 
C € [0, 1) when t — * 00. We note that Q = k = 
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FIG. 6: The feasible regions for {(, k) if the MSD is 
diffusive at small time. 



is the special case of the Dirac delta function kernel. 
The fractional generalized Langevin equation considered 
in Section IIII Bl corresponds to k = (, whereas the ex- 
tended fractional generalized Langevin equation consid- 
ered in Section IIII CI corresponds to K = 0. For the 
Laplace transform of the dissipative memory kernel j(s), 
one finds that j(s) ~ when s — > and 'j(s) ~ s~ K 
when s — > oo. The solution x(t) can be written as (|28[). 
where the small-i and large-i asymptotic behaviors of 
the functions K\{t) and K 2 (t) are still given by ([50]) 
and (HIT]) . The generalized fluctuation-dissipation theo- 
rem (|32|) then implies that similar to j(t), Cp{t) ~ £ K 
when i — > and Cf(£) ~ when t — > oo. One can 

then deduce that when t — > 0, the MSD (|29|) behaves like 



A 2 W 



if 2a + K> 1, 




FIG. 7: The feasible region for (£ «) if the MSD is ballistic 
at small time. 



or (£, k) satisfies (see FIG. [7|) 

- 2( > -1/2, k < 1, 



1 ^ 3 
2 <C< 4' 



(37) 



and 



<* = --<;, (3 = i. 

Finally we would like to comment that if the dissipative 
memory kernel is a finite sum of power-law functions, 
then k < Ci an d the situations (f3"6")l and (f3"T)) would not 
appear. 



if 2a + K < 1. 



IV. 



CONCLUDING REMARKS 



The large-time asymptotic of the MSD is still given by 
(fM)) . After some analysis, we find that the short and long 
time properties of SFD are satisfied if (C, k) satisfies (see 
the first graph in FIG. [5]) 

C>0, K -2C>i, k<1, 



and the values of a and (3 are 
1 



1 



- 4 C ' "=2' 
or (C, k) satisfies (see the second graph in FIG.[H1) 

C>0, k>0, K -2C<i, k + 2C<^, 
and the values of a and f3 are 



(36) 



4 2' P 8 



On the other hand, the MSD is ' 
t if KE [0,1), (3 = 1, C = 1/2 and 



4 2 
t 2 (ballistic) at small 



, 1 1 - « 

a > max < — , 

^ 4' 2 



We have proposed some stochastic processes that can 
describe the basic properties of SFD. In the first case, 
we find that the replacement of the Hurst exponent in 
fractional Brownian motion by a two-valued step func- 
tion can be used to model SFD that is normal diffusive 
at small time and sub-diffusive at large time. Since frac- 
tional Brownian motion can be considered as a solution of 
a special case of fractional Langevin equation, it is nat- 
ural to consider modeling SFD by fractional Langevin 
equation. We have discussed in detail the cases where 
the dissipative memory kernel is a Dirac delta function, 
a power-law function as well as the combination of these 
two. For each of these cases, we find that there is a 
range of the parameters where the corresponding frac- 
tional Langevin equation can be used to model process 
whose MSD behaves like ~ t when t is small and behaves 
like ~ \ft when t is large. This range has been determined 
explicitly. The corresponding range of the parameters for 
which the MSD of the process behaves like ~ t 2 instead 
of ~ t when t is small is also determined. 

One weakness of usig SFBM to model SFD is that the 
change of scaling exponent is not continuous so that the 
transition from the normal diffusion to the anomalous 
subdiffusion occurs abruptly. This weakness is overcome 
by various fractional Langevin models. However, we only 
use the short and long time behaviors of the stochastic 
processes to determine the model. For a more complete 
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description of SFD, it is necessary to take into account 
the boundary conditions at the ends of the single-file sys- 
tem, as these would become relevant if the system is of 
finite extension. Furthermore, one may also have to con- 
sider the interaction between the diffusion particles with 
the wall of the single-file system. Thus, the evolution 
of the particles in the single-file system during the in- 
termediate time interval becomes relevant. We hope to 
incorporate the physical mechanism of SFD in order to 
obtain a more realistic model of SFD in our future work. 
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APPENDIX A: THE FUNCTIONS ifi(t), K 2 (t) 
AND THEIR ASYMPTOTIC PROPERTIES 



In this appendix, we give the details of the analysis 
of the asymptotic properties of the functions K\ (t) and 
K 2 (t) defined as the inverse Laplace transform of 



where Ko(t) is the inverse Laplace transform of 

1 



1 + Xs- a + xs~ c ~ a ' 
which has an integral representation 



-1 pioo 

Ko(t) — / - 



ds 



e ist 

2tt J 1 + \(is)- a +x(is)-<~ a 



1 



ds 



2tt J 1 + \(-is)- a + x(-is)- c ~ a 



ds 



ds 



Im 

" 7r J 1 + \e™ a s- a + xe irr «+ a ^s^i- a 

Im r°° s^ +a e- st 

ds 



Now for the large-t asymptotic behavior of Ko(t), we have 



(f y + Ae"« (f y + X e i7r (C+«) 



-ds 



and 



s a+0 + X s -a + xs -C-a) s 0+l (1 _|_ \ s -a _|_ ^g-f- 

By taking Laplace transforms, it is easy to check that 



Im 



k=0 



3=0 

ft 



- J x J 



r (ah Tcf+a+Jj ' 

jak+Cj + 8 

Therefore, as t — > 0, 

Ki(t) ~ t a+fi ~ x and K 2 (t) ~ t (i . 

This behavior is independent of the coefficients A and x- 
For the large-i asymptotic behaviors, notice that 



Ki{t) 



K 2 (t) 



1 



r(a + 13) J 



1 



(t-u) a+fi - 1 K (u)du, 
(t - uf K (u)du, 



(Al) 



t -Ck- a j 8 <;k+ a j+c+<x e -s ds 



fc=0 



3=0 



x r(Cfe + aj + c + a + i)r Cfc - ai 

This implies that generically, 

K (t) ~ r 1 -^- , t -► oo. 
We then obtain from ()A1[) that as i — > oo, 
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